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Abstract 

We study color superconductivity in QCD at asymptotically large chemi- 



\o 
o . 

Q>^ , cal potential. In this limit, pairing is dominated by perturbative one-gluon 

exchange. We derive the Eliashberg equation for the pairing gap and solve 
, this equation numerically. Taking into account both magnetic and electric 

(— I i gluon exchanges, we find A ~ g~^ ex]3{—c/g) with c = 37r^/\/2, verifying 



a recent result by Son. For chemical potentials that are of physical interest, 
/x < 1 GeV, the calculation ceases to be reliable quantitatively, but our results 
suggest that the gap can be as large as 100 MeV. 
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1. The behavior of matter at very high baryon density but small temperature is of interest 
in connection with the physics of neutron stars and heavy ion collisions in the baryon rich 
regime. Moreover, it has been realized that matter at very high density exhibits many non- 
perturbative phenomena, such as a mass gap and chiral symmetry breaking, in a regime 
where the coupling is weak and systematic calculations are possible. 

At very high density the natural starting point is a Fermi sphere of quarks. The cor- 
responding low energy excitations are quasiparticles and holes in the vicinity of the Fermi 
surface. Since the Fermi momentum is large, asymptotic freedom implies that the interaction 
between quasiparticles is weak. However, as we know from the theory of superconductiv- 
ity, the Fermi surface is unstable in the presence of even an arbitrarily weak attractive 
interaction. In QCD, the attraction is provided by one-gluon exchange between quarks in 
a color anti-symmetric 3 state. QCD at high density is therefore expected to be a color 
superconductor [||,|1. 

A particularly interesting case is QCD with three flavors. In this case the most favorable 
type of pairing involves the coupling of color and flavor degrees of freedom, color-flavor- 
locking [^ . This implies, among other things, that all gluons acquire a mass and that chiral 
symmetry is broken. We have argued that in the color-flavor-locked phase not only universal 
features, in particular the symmetry breaking pattern, but also many non-universal proper- 
ties, such as the spectrum of low-lying states, exactly match the expectations for hadronic 
matter at low baryon density [^]. This means that nuclear matter at low density might be 
continuously connected to quark matter at high density, without any phase transition. 

It also means that at very high density many interesting properties of hadronic matter, 
such as the magnitude of the chiral condensate, can be calculated in weak coupling perturba- 
tion theory. The first step in such a program is the calculation of the superconducting gap. 
This calculation was first attempted by Bailin and Love 0. They used a schematic IR cutoff 
in the gluon propagator. In their treatment, the gap is of the form A ~ fiexp{—c/g^), where 
c depends logarithmically on the IR cutoff. Recently, the problem was revisited by Son [Q, 
who argued that the gap is dominated by magnetic gluon exchanges, and that the infrared 



behavior is regulated by dynamic screening. He obtained A ~ fig~^ exp^—STr"^ / {\^g)) . 

Our purpose in the present work is to rederive and strengthen this result, and to deter- 
mine the overall numerical coefficient. 

2. In order to derive a gap equation, we follow the standard Nambu-Gorkov formalism 
and introduce a two component field \E' = {ip,ip'^). The inverse quark propagator takes the 
form 

I 'A+ d-m A » 
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where A = 7oA'''7o. The gap is a matrix in color, isospin, and Dirac space. In the following 
we consider the case of two massless flavors. We will assume that the gap is anti-symmetric 
in both flavor and color, and has total angular momentum zero. In the case of short range 
interactions, this assumption can be justified from our study of the renormalization group 
equations for a general four-fermion interaction 0,H. The one-gluon exchange interaction is 
long range, and other forms of pairing might take place. In particular, since the interaction 
is dominated by almost collinear scattering, one might expect higher partial waves to play 
a role 0. In the following, we will concentrate on the total angular momentum zero gap. 

We also assume that the gap has positive parity. One-gluon exchange does not distinguish 
between pairs of positive or negative parity . This degeneracy is lifted by instantons, which 
favor the positive parity channel |p!0| , pn] ,|7| . At large chemical potential instanton effects are 



exponentially suppressed. In the following, we will therefore assume that the only instanton 
effect is to determine the parity of the gap. 

In addition to that, we neglect quark mass effects and chiral symmetry breaking LR 
condensates. As shown in there is no BCS instability in the case of pairing between left 
and right handed quarks. The formation of LR condensates is therefore suppressed by m/fi. 
The form of the gap matrix is then |5|,p!2 

A^^g) = i\2nr,),,C^, (Ai(go)^(l + a ■ q) + A^iqo)^^ - « ■ g)) • (2) 

In the weak coupling limit, we can replace a ■ q hj the unit matrix using the equations 



of motion. In this limit, only Ai survives. We will see this explicitly from the solution 
of the gap equation presented below. We have neglected the dependence of the gap on 
the magnitude of the momentum, but kept the dependence on frequency. The dependence 
on momentum can be dropped because, in the weak coupling limit, all momenta are close 
to the Fermi surface. For short range interactions, the dependence on frequency can also 
be neglected. This is not the case here. Because long range interactions are important, 
retardation effects cannot be neglected. 

The self energy in the Nambu-Gorkov formalism obeys the Dyson-Schwinger equation 

i 

S(A:) = -^g' J -0flS{q)TlDfM - k). (3) 

Here, S(A;) = —{S^^{k) — SQ^{k)) is the proper self energy, F^ is the quark-gluon vertex and 
D°^^{q — k) is the gluon propagator. To leading order in the perturbative expansion, we can 
use the free vertex 



■pa 



/^7^AV2 ^ 



(4) 



\^ -(7,AV2)^; 

We will study the importance of vertex corrections below. To leading order, we can also 
neglect the diagonal part of the proper self energy, that is the fermion wave function renor- 
malization 0. In this case, @ reduces to an equation for the gap matrix, 

A(fc)=^//^(7.y) S,,{q)(^^^^^D,^{q-k). (5) 

Here, S2i{q) is the 21-component of the fermion propagator in the Gorkov representation. 
S2i{q) is determined from the inverse of (|l|). We have 

Inserting the ansatz (|^) for the gap gives 

„,> 1,, „ J Ai(l-a.j) , A,(l + S-q} \ 

&(,) = -5(A.r,C-,,) I , _ 11^^ _ ^j, _ ^, + ^, _ ^1^^ ^ ^j, _ ^J . (7) 



Both the RHS and the LHS of the gap equation are proportional to T2, so the flavor structure 
simply drops out. The color coeflicient is given by 
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where we have used the Fierz identity (A")jj(A")fe/ = —2/NcSij6ki + 26ii6jk and the factor 1/4 
conies from the color generators f^ = X°'/2. Projecting (^ on Ai^2 gives two coupled gap 
equations 



2ig^ f d^q (1 / ^ LaA ^i(go) 
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where the two signs oi a ■ k on the RHS correspond to Ai and A2 on the LHS. 

We now must specify the gluon propagator. The gluon propagator in a general covariant 
gauge is given by 

P^ P^ a a 

D,M = -T^ + -r^ - e^ (10) 

q^ — G q'^ — F q* 

where D and F are functions of go and \q\ and the projectors -P^^ are deflned by 

Pl,=5,,-m,. P^o = P^. = o, (11) 

It contains the gauge parameter ^, which must not appear in physical results. In the weak 
coupling limit, go is small as compared to |g|. In this case we can expand the projectors 
P^j^ ~ 5^q5i,q and qtqj/q^ — (jiqj- The gap equation now becomes 

2tg^ f dS \ Ai(go) ( I'l^-q , \ + \k-q 
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There is a similar equation for A2 in which the two terms in the round brackets are inter- 
changed. Only the flrst term term in ([l3|) has a singularity on the Fermi surface. In the 



weak coupling limit, we can therefore drop the second term, and we are left with an equation 
for A{pq) = Ai(po)- This equation is independent of the gauge parameter ^. The second 
gap parameter A2 is not suppressed in magnitude. However, A2 does not lead to a gap on 
the Fermi surface, and its value is gauge dependent. 

We should note that the fact that the gap is gauge independent in the present weak- 
coupling approximation is a consequence of the fact that the gap is determined by the 
scattering of quarks that are almost on shell. For on-shell quarks, the fact that the gauge 
dependent part of the propagator does not contribute follows directly from the equations of 
motion for the quark fields. 

For large chemical potential the integral over q is dominated by momenta in the vicinity 
of the Fermi surface, |g| ^ /i and go *^ f^- We can expand all momenta as q = qp + 1, where 
qp is on the Fermi surface, and / is orthogonal to it. Asymptotically, |/| ^ {qpl and the 
integration measure becomes dqo jj? did cos dcj). We also have \q — k\ ~ v2Ai(l — cos^). The 
integral over is performed trivially. We analytically continue to imaginary go, and perform 
the integral over / by picking up the pole in the diquark propagator. We find 

Q^ f f ( - — - COS Q 

i + |cos^ \ A(go) 



1 _ cos ^ + (F + (po - go)')/(2/i2) ; y^g2 + A(go)2' 

The integral over cos^ is dominated by small ^, corresponding to almost collinear scattering. 
It is therefore important to take medium modifications of the gluon propagator at small 
momenta into account. For go ^ g ^ and to leading order in perturbation theory we have 

F = 2m^ G = |m^^, (15) 

with 171^ = Nfg'^fi'^/{An'^). In the longitudinal part, m|, = 2m^ is the familiar Debye 
screening mass. In the transverse part, there is no screening of static modes, but nonstatic 
modes are modes are dynamically screened due to Landau damping. In our case, typical 
frequencies are on the order of the gap, go — A. This means that the electric part of 
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the interaction is screened at qE — rrij^ whereas the magnetic interaction is screened at 
qm ~ (7r/4 ■ m|)A)^/3_ 

Asymptotically, qm -C Qe, and magnetic gluon exchange dominates over electric gluon 
exchange. We therefore begin by analyzing the gap equation taking into account the mag- 
netic part of the interaction only. We will also approximate cos ^ ~ 1 in the denominator 
and drop (go — Po)^ in the denominator. All of these terms will be reinstated later. The 
integration over cos^ is now straightforward. We have 

A(Po) -^Jd,, log f 1 + ^^^ -^^- (16) 

If we are only interested in the leading exponential behavior of the gap we can drop the 
numerical factors and the powers of g in the logarithm. We then arrive at 

A(Po) = 4-J'^^ ^°S f HS) T#^^' (17) 



which is the equation discussed in the appendix of Son's paper 0. This equation was 
derived from the on-shell quark-quark scattering amplitude. What we have shown here is 
that one can indeed derive this equation from the Dyson- Schwinger equation in the weak 
coupling limit, and that the result is independent of the gauge parameter. Son also derives 
an approximate solution to this integral equation, 

A,pp.(po)=Aosin(^-^log^^jj, Po>Ao, (18) 

with Aq = /iexp(— 37r^/(-\/2(7)). The approximations involved are expected to reproduce the 
correct coefficient in the exponent, but do not fix the prefactor. 

3. We have therefore solved the Eliashberg equation ([I7|) numerically for different chemi- 
cal potentials. We have used the one-loop running coupling constant evaluated at the Fermi 
momentum pp = yU- This is an average over the momenta of the exchanged gluons, which 
are in the range [qmi 2yu]. Without a higher order calculation one cannot fix the scale in the 
running coupling. We will see that the preexponential factor in the final result behaves as 
g~^ . This factor is almost optimal to give a remarkably weak scale dependence. 



The result for the function A(po) for /i = 400 MeV and fi = lO^'' MeV is shown in Fig. |l|. 
The sohd hne is the numerical result while the dashed line shows the approximate solution 
(0), rescaled by an overall factor c, A(po) = cAapp. (po/c). At yU = 10^° MeV, g ^ 0.67 and 
Son's solution is in excellent agreement with the exact result, up to an overall factor c ~ 2. 
At /i = 400 MeV the coupling is significantly bigger than 1, (^ — 3.43, but the approximate 
solution is still qualitatively correct. 

The scaling of the maximum gap with the chemical potential is shown in Fig. ^. The 
solid line is the numerical result and the dashed lines correspond to c(7~'^/i exp(—37r^/(-\/25f)) 
with k = 0, . . . , 5. We observe that the k = curve provides an excellent fit to the data 
even for small chemical potentials. Again, the overall coefficient is c ~ 2. 

Let us make a few observations at this point. First, we note that the use of perturbation 
theory to determine the dynamic screening is self consistent. Since A ~ /xexp(— const. /^f), 
the gap grows as /i — i> oo and qm ^ ^qcd- Second, we note that it is essential to keep 
the frequency dependence of the gap. For small frequencies A(po) varies over scales on the 
order of po ~ Aq itself. Therefore, A{pq) cannot be replaced by a constant. Were we to 
approximate A{pq) ~ Aq, as in W^, we would obtain a gap equation for Aq that has the 



correct double logarithmic structure and gives Aq — /iexp(— const. /(?), but the constant in 
the exponent would not be correct. 

Finally, we note that it is easy to see what taking into account the numerical coefficients 
and the factor g^ in equation ([T6|) will do. Any numerical factor inside the logarithm can 
be absorbed by rescaling the frequencies. Therefore, if Aapp.(Po) in (0) is an approximate 
solution to (0), then A'(po) = cAapp.(po/c) with c = GAn/^Nfg'^) is an approximate solution 
to (p!6D. This can also be seen from Figs. ^ and | where we show the numerical solution to 



the Eliashberg equation (p!6| ) for the superconducting gap from magnetic gluon exchanges. 
Asymptotically, the solution is well described by the function A'(po) with c ^ 175g~'^. 

We now come to the role of electric gluon exchanges. We include the second term in ( p!^ 
with F = rn^). We again use the approximation cos 9 ~ 1 in the numerator and drop the 
(go — PqY term in the denominator. Let us note that in the forward direction, electric and 
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magnetic gluon exchanges have the same overall factor. Performing the integral over cos 9, 
we find 



/qi + Mqof 

where the factor 3/2 in front of the second term comes from the difference between dynamic 
screening, q^ ~ |^^^'^, and static screening, qE ^ \q\- In the weak coupling limit we again 
expect the solution to be of the form cAapp(po/c) with 

c = 102AV2tt^NJ^^^ g-^ = 256TT^g-^ - 2.5 ■ lO^g'^ {Nf = 2). (20) 



We can compare this prediction to our numerical results, obtained from solving (|14[). In 
this equation, we take into account both electric and magnetic gluon exchanges. We also 
keep the cos 6' dependence in the numerator, and the terms (go — Po)"^ in the denominator. 
Finally, we use the exact form of G and F in the hard dense loop approximation, 

F = 2,„<fl-^Q,felY Q„(.) = i,ogfi±i-). (21) 



Q^ \ 1^ \\<l\/J' 2 Vx — 1 



Cr = m -— 

\q\ 



2\ 



IMM 



(22) 



This takes into account that there is no dynamic screening for |g| < go- The numerical 
results are shown in Figs. || and |. Asymptotically, the gap is well described by cAapp(po/c) 
with c ~ 1.4 ■ 10^g~^. We notice that for /i = 10^'^ MeV the solution has a 'knee' at po — 10^ 



MeV. This comes from the fact that for frequencies po > JNf/{87r)gn the retardation terms 
~ {pq — go)^ dominate over screening. In this regime, the solution is of the same form, but 
the scale factor is different. 

Overall, the scaling with 5f~^exp(— 37r^/(v^5f)) is clearly visible, though not quite as 
impressively as in the case with magnetic gluon exchange only. For chemical potentials that 
are of physical interest, /i < 1000 MeV, the gap reaches Ao — 100 MeV. We should caution, 
however, that in this regime g ~ (2—4), and higher order corrections are probably important. 
Nevertheless, it is gratifying to see that the order of magnitude of the result agrees with 
previous calculations 0,Q based on more phenomenological effective interactions, which 
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were normalized to the strength of chiral symmetry breaking at zero density, rather than 
the calculable asymptotics of the running coupling. 

4. There are a number of questions that will need to be addressed in a more complete 
calculation. First, we have concentrated on the case Nf = 2. For Nf = 3, there are two order 
parameters, corresponding to the color antisymmetric and color symmetric components of 
the color-flavor locked state. This is only a minor complication, since there is only one 
combination that survives in the weak coupling limit. 

A more complicated issue is the role of the Meissner effect. For Nf = 2, the dominant 
order parameter only breaks color SU{3) -^ SU{2), and all gluons that contribute to pairing, 
except for one, live in the unbroken part of the gauge group. In the case of Nf = 3, the Higgs 
mechanism is complete and all gluons acquire a mass. At zero momentum and frequency, the 
screening mass is on the order of m'^ ~ g^fi^, much larger than the dynamic screening scale 
qm- At finite momentum transfer, on the other hand, the screening mass is m'^ ~ g'^n'^A/\q\ 



Ij], which is of the same form as the dynamic screening effect. The Meissner effect will 
therefore not affect the dependence of the gap on the coupling constant, but it will affect 
the numerical coefficient. 

Finally, one has to address higher order corrections to the perturbative result. In par- 
ticular, one would like to know what the functional form of the corrections is, and whether 
the applicability of perturbation theory requires g < 1, or some weaker condition like g < ir. 
We have already mentioned wave function renormalization as one source of higher order 
correction [0. Another issue is vertex corrections. The vertex correction generated by hard 



dense loops is 15 






r;{pi,P2) =9-[l,-m)l -7-^6-^ 1 ' (23) 



where K = {i, k) is a light like vector and mi = 5f^/i^/(67r^). We can insert this correction 
into the gap equation (j^). We find that the coefficient of the magnetic gluon exchange is 
modified as 
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3 _ icos^) ^ (- - -cos^) + ^ + -(1 -cos^K f^ ^ — (24) 

In the forward direction cos 6' ~ 1, which dominates the gap equations, this is just a higher 
order correction. Vertex corrections have a (fi/qo) enhancement in the backward direction 
cos^ ^ —1, but the integral over cos^ is finite as go ~^ 0. Vertex corrections in the magnetic 
part therefore do not modify the asymptotic form of the gap. The same is true for vertex 
corrections in the electric part of the interaction. 

5. In summary, we have performed a perturbative calculation of the superconduct- 
ing gap in two fiavor QCD at very high density. We find that the gap scales as Aq — 
2567r^/i5'^^exp(— 37r2/(v^5f)), where the overall coefficient is correct up to a factor of order 
one. In the physically interesting regime n < 1 GeV, the gap is on the order of 100 MeV, in 
agreement with earlier calculations based on instantons or schematic interactions adjusted 
to the size of the chiral condensate at zero density. 

Acknowledgments: We would like to thank D. Son and K. Rajagopal for useful discus- 
sions. 
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FIGURES 
FIG. 1. Solution of the Eliashberg equation (|l^) as a function of imaginary frequency qq. The 

upper and lower panels show the solution for ^u = 400 MeV and fi = 10^*^ MeV, respectively. The 

solid lines show the numerical solution and the dashed lines shows the approximate solution (p!8|), 

scaled to the same value of the gap. 

FIG. 2. Dependence of the gap on the chemical potential for the solution of the Eliashberg 
equation (|l3)- Here, ^(/u) is taken to run according to the one-loop beta function. The dotted 
curves show the functions g^^ exp(— (37r^)/(\/25)) for fc = (top), . . . , 5 scaled to the value of the 
gap at the maximum chemical potential. 

FIG. 3. Same as figure |l| for the solution of the Eliashberg equation with magnetic gluon 
exchange only, see equ. (Wt 



FIG. 4. Same as figure |2| for the solution of the Eliashberg equation with magnetic gluon 
exchange only. 

FIG. 5. Same as figure || for the solution of the Eliashberg equation with magnetic and electric 
gluon exchanges, see equ. ([l^. 



FIG. 6. Same as figure g for the solution of the Eliashberg equation with magnetic and electric 
gluon exchanges. 



13 




10 

q„ [VeV] 





- ' 


' ' 


■ 


0.0008 


~~~*"--«<^ 




- 


0.0006 


- 


\^ 


- 


^.0004 




\ ■ 
\ ■, 


- 


0.0002 


- 


\^ 













0.00010.001 0.01 0.1 1 10 100 1000 

q„ |MeV] 

FIG. 1. 



lO* 10= 10" 10' 1D» 10" 10'° 



14 



10 



> 1 



0.01 



0.001 



" 


1 1 1 IMIj 1 1 1 MMIj 1 1 M MM 1 1 1 


nil 1 1 1 1 1 III 


1 1 1 1 nil 1 1 1 1 1 III 1 1 1 1 1 III 


1 1 1 1 1 III 1 1 1 1 1 II 


1 


- 


\ '.1 
\ '.\ 
\ 'I 
\ \ 

\ A 

\ 'A 

^ ''\ 
\ \ \ 

\ ^ \ 








- 


' 


\ \ \ 








" 


' 


\ V \ 








" 


- 


\ \ ^\ 








- 


- 


\ \ X \ 

\ \ \ \ 
\ \ \ \ 








- 


- 


\ \ \ -^ 








- 


- 


\ \ \ 'X 








- 


r 


^ \ \ \ 'X 
\ \ -. "■ 

\ \ V, 

1 ^ 1 "- 1 "^ 

1 Mil -u III 


iiiil^ r-i oaiiiil 


_--^ 


1 




1 1 iiiiJ- -r iTi ij^iX^ i>- 1 



1000 



104 105 



10^^ 10'' 108 10" 101° iQii 101 

/i[McV] 



FIG. 2. 



15 



10 



5 - 



10 



q„ iMeVj 



100 



0.12 



0.1 



0.08 



0.06 - 



0.04 - 



0.02 - 




16 



10 



> 

CD 



0.1 



- 


I 1 M[ 1 1 1 Mill 1 1 1 1 Mil 1 1 1 1 nil 1 III 

1 \ 

\ ' * 
1 1 \ 

'',\ 1 \ 
'\ I \ 


1 1 1 1 1 nil 


1 1 1 1 III 1 1 1 1 1 III 1 1 1 1 1 III 1 1 1 1 1 1 


v^ 


- 


A \ ^ 
\ 'A \ \ 
\ 'A \ \ 






- 




\ A \ \ 
\ A \ "^ 

\ \ \ \ 








7 


\ \ \ ^ 






7 


- 


\ A "^ 






- 


- 


\ \ \ \ X 




.^j/ 


- 


- 


\ \ \ \ X 

\ \ \ ^ "" 






- 




\ \ \ "-^ 


~ ^ _ 








\ \ V -^ 




^ -^ ^^^ ^/ 






\ \ ^-\ 


^ — — 


C^;/"'' 




: 


\ ^ ^ ^'^'^ 


.,„^^,— -— " 


: 


- 


\ 1 "x 

i-t-i mil. 1 1, 1 1 1 iTtll.ll 


1 1 1 ^inrl 


\A V 





1000 



104 105 IQS 



10'' 108 10" 101° iQii 101 

/j,[MeV] 



FIG. 4. 



17 



100 



80 - 



60 



40 



20 



100 



q„ |MeV| 




100 1000 



10" 10'' 10= 

q„ [MeV] 



FIG. 5. 



100 



> 

(1) 



10 




iiJ I I I I 



1000 



10* IC 106 



10^ IQB 

^t[MeV] 



109 iQi" 10" 101 



FIG. 6. 



19 



